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We study the phase structure of a dilute two-component Fermi system with attractive interactions 
as a function of the coupling and the polarization or number difference between the two components. 
In weak coupling, a finite number asymmetry results in phase separation. A mixed phase containing 
symmetric superfluid matter and an asymmetric normal phase is favored. With increasing coupling 
strength, we show that the stress on the superfluid phase to accommodate a number asymme- 
try increases. Near the infinite-scattering length limit, we calculate the single-particle excitation 
spectrum and the ground-state energy at various polarizations. A picture of weakly-interacting 
quasi-particles emerges for modest polarizations. In this regime near infinite scattering length, and 
for modest polarizations, a homogeneous phase with a finite population of excited quasi-particle 
states characterized by a gapless spectrum should be favored over the phase separated state. These 
states may be realized in cold atom experiments. 



Recent experiments on cooled Fermi atoms [B S S Q 
and theoretical developments in dense QCD [f| have mo- 
tivated renewed theoretical interest in Fermion superfiu- 
ids. In Fermi systems it is well known that attractive 
interactions destabilize the Fermi surface. This insta- 
bility is resolved by the BCS mechanism characterized 
by pairing between spin-up and spin-down particles with 
opposite momenta on the Fermi surface. The system 
exhibits superfluid properties and has a gap in the ex- 
citation spectrum. Measurements have been performed 
probing the equation of state (EoS) |2j and the pairing 
gap m the strongly-interacting regime. Quantum 
Monte Carlo (QMC) methods have also been employed 
to examine these properties 0, S @ • 

These studies have primarily addressed the unpolar- 
ized system. In contrast, the ground state properties of 
an asymmetric two-component system remain unclear as 
several competing states have been proposed. These in- 
clude: a gapless superfluid [ToL HH ll^. fili. Hif : a mixed 
phase consisting of BCS and normal components |l5j| : 
and the Larkin, Ovchinnikov, Ferrel and Fulde (LOFF) 
phase where the order parameter acquires a spatial vari- 
ation [lfl Il7| which in three dimensions may manifest 
as crystalline state 01- The competition between these 
states is especially relevant to understanding the phase 
structure of dense quark matter JjJ. In this letter we 
address several issues relating to the ground state of an 
asymmetric two-component Fermi system in the strong 
coupling regime. Using both mean-field theory and QMC 
we find that in strong coupling and at finite polarization 
the stress on the BCS state to accommodate the number 
asymmetry increases. In the intriguing strong-coupling 
regime near /cpa = oo, a superfluid state with non-trivial 
gapless excitations may be favored. We also find that 
these quasi-particles are weakly interacting. 

Two- component Fermi System: We consider a sys- 
tem consisting of non-relativistic spin-up and spin-down 
Fermions at finite polarization. A short-range potential 



(range -C inter-particle distance) between spin-up and 
spin-down particles characterizes the interaction. The 
interaction between spin-up particles (or spin-down par- 
ticles) is considered to be negligible compared to the in- 
teraction between up and down spins. We study both 
weak coupling and strong coupling limits to explore the 
phase structure from the BCS-like regime characterized 
by pairing at the Fermi surface to the Bose Einstein con- 
densate regime characterized by bound bosonic states. 
The Hamiltonian (Grand-Canonical) is given by: 
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where g is the effective four-fermion interaction whose 
strength at low energy is determined by the two-body 
scattering length a. For the two component system, the 
spin-up and and spin-down chemical potential may be 
written as ^ = /i+8fi and fj,^ = /i—Sfi, respectively. The 
density n — n-\ + n± determines fi and the polarization 
density Sn — — nj, determines Sfi. In trapped atom 
experiments this regime can be accessed by adjusting the 
population of the two species. 

Initially, we consider a Fermi system in the "univer- 
sal" regime characterized by an infinite scattering length. 
When a — oo , the interaction does not present a dimen- 
sionful scale. Consequently, the energy density, pressure 
and chemical potential of the unpolarizcd system are re- 
lated to that of the Fermi gas by the relation e(a = oo) = 
£ €fg, P( a = oo) = £ Pfg, an d A* = £ k F /2m, respec- 
tively. Quantum Monte Carlo studies in small systems 
with 12-20 particles have determined this numerical co- 
efficient to be £ = 0.44(2) 0. The system also exhibits 
a gap in the excitation spectrum, numerical studies indi- 
cate that the gap A = 0.95(5)E FG = fin with (3 s» 1.4. 
Below we present new calculations for larger system sizes, 
which are useful to more fully explore the dispersion of 
the single-particle excitations. For these larger systems 
we find £ = 0.42(1) and A = 0.84(5)E FG . We note that 
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the value for £ is in agreement with experimental studies 
described in Ref. 0- 

Our primary interest here is the spin-polarized sys- 
tem. Earlier work by Bedaque, Caldas and Rupak based 
on BCS mean-field theory showed that a finite polariza- 
tion would lead to phase separation |15|, A heteroge- 
neous mixed phase consisting of an unpolarized super- 
fluid state coexisting with a partially polarized Fermi 
gas state was shown to have lower Free energy than 
the homogeneous gapless superfluid phase (also called 
the breached-pair phase) suggested earlier by Liu and 
Wilczek 0| . More recently, Forbes, Gubankova, Liu, 
and Wilczek 0] have shown that the this phase may 
be stabilized by finite range interactions and different 
masses for the two species. Here, though, we consider 
the equal mass case with short-range interactions, exam- 
ining the phase structure at stronger coupling. 

Normal- Superfluid Mixed Phase: Polarizing the super- 
fluid state is disfavored due to the presence of a gap in 
its excitation spectrum A heterogeneous state contain- 
ing normal and superfluid phases provides an alternate 
route to accommodate a finite polarization. Here the ex- 
cess spin-up particles could reside in the normal phase. 

Phase coexistence is possible between states separated 
by a first order transition if at fixed chemical potential 
they can have the same pressure. Hence we require 



Substituting the QMC results for £ and f3, we find that 
Sfx/A = 1.00(5). This suggests that the BCS state is 
close to the polarization threshold; and that by tun- 
ing the coupling strength cold-atom experiments in traps 
containing two-Fcrmion species should be able to explore 
the large Sfi/ A regime. Further, we find that Sfi > fi 
- indicating that the polarization in the normal phase 
is maximal. It is intriguing that we find A ~ Sfi at 
kpa — oo. 
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FIG. 1: Variation of Sfi/ A with increasing coupling. 



Pressure equilibrium uniquely determines Sfi for a given 
fi. We emphasize that at fixed fi, Sfi in the mixed phase 
does not change with polarization - it is driven to lie 
exactly at the first order transition point satisfying Eq.^ 
An increase in net polarization is accommodated by an 
increase in the volume fraction of the normal phase. 

We observe that Sfi is the energy required to intro- 
duce a spin-up particle into the normal component of 
the mixed phase and A is the corresponding energy in 
the superfluid component (where the energy is measured 
with respect to the chemical potential fi). In the weak 
coupling BCS regime, the mixed phase is characterized 
by Sfi/ A = 1/ V2 and the BCS state remains unpolarized. 
When Sfi > A the BCS state will acquire a finite polar- 
ization resulting in gapless excitations in its spectrum. 
Otherwise^a finite polarization will result in phase sepa- 
ration 15,20]. Thus, at infinitesimal polarization, <S fi— A 
is the energy difference per unit polarization between the 
homogeneous gapless phase and the mixed state. These 
observations motivate us to examine how the ratio Sfi/ A 
changes with coupling strength. 

We first analyze the situation at kpa — oo. Using the 
EoS results of Ref. Q and conditions of Gibbs equilib- 
rium (Eq. ^) we find that the ratio 
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The gap increases with coupling but Sfi increases at a 
faster rate. The increase in Sfi is driven by the pressure 
equilibrium condition. This trend can be demonstrated 
by using mean field theory in the weak coupling regime 
where A(kpa) = 8/xexp (— 2 + (7r/2\/2fcFa)) • Earlier 
studies of the mixed phase ignored the leading order 
kFa corrections to the pressure and chemical poten- 
tial • However, including these corrections is straight- 
forward [il| . We find the pressure in the normal and BCS 
state is given by 
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respectively, where Pa = mkpA 2 /(4-7T 2 ) is the pairing 
contribution to the pressure. The chemical potentials 
also receive corrections at order k-pa. In the normal phase 
MTU) = fc FT(i) /2m+ (iira/m) fc Fm) /(67T 2 ). 

Gibbs equilibrium condition determines the variation 
of Sfi as a function of 1/kpa at fixed fi. While we expect 
the mean field result to be valid only for small coupling, 
the results which are shown in Fig. show the expected 
trend - an increase in Sfi/ A with increasing coupling. The 
QMC result at the universal point k-pa = oo confirms this 
increase. The behavior in the extreme BEC limit is also 
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shown. In this region, we may use mean field theory to 
calculate the self energy of a spin-up Fermion in the BEC 
phase. This is given by Abec = ^naBF^/m, where 
cibf — l-2a is the Fermion-Boson scattering length [2^ . 
ub is the Boson density and fh is the reduced mass of 
the Fermion-Boson system. The Boson-Boson scatte ring 
length is also known and is given by obb — 0.6a |22fl . 
We use this to calculate the pressure of the BEC at lead- 
ing order in the n^aBB expansion 21]. Pressure and 
chemical equilibrium then uniquely determine the chem- 
ical potential of spin-up Fermions in the Fermi gas phase. 
In agreement with earlier work by Viverit, Pethick and 
Smith [2^ we find that the homogeneous BEC phase 
easily accommodates a finite polarization in the dilute 
regime - as evidenced by 5/x/Abec 3> 1. 
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FIG. 2: The quasi-particle spectrum above the superfluid 
phase at tea = 00. For reference, the quasi-particle spec- 
trum in BCS theory for k^a — 00 is also shown. 

Quasiparticle Dispersion and Polarized Superfluid 
State: In addition to calculating the paired ground-state 
energy and the superfluid gap, we have performed QMC 
calculations of the homogeneous superfluid phase to ex- 
amine the quasiparticle dispersion as a function of mo- 
mentum. In addition, we have calculated the ground- 
state energy of systems at finite polarization to exam- 
ine the interaction between these quasi-particles and de- 
termine if the gapless phase can support a macroscopic 
polarization. The methods used are identical to those 
employed earlier and the same (finite-range) cosh poten- 
tial was used [1,0- As before, we expect the small but 
finite-range of the potentail to have a small effect upon 
the ground-state energy and the gap. Here we employ 
somewhat larger system sizes, however, ranging from 54- 
66 particles rather than the earlier studies from 12 to 20 
particles. The larger system sizes allow for a somewhat 
finer momentum grid which is useful for examining the 
dispersion. 

For the dispersion calculations, we place an unpaired 



spin in a state of definite momenta k = (n x ,n y , n z )2ir/L, 
where L is the (cubic) box length and the ni are integers 
describing the momenta in each coordinate. The BCS 
plus unpaired particle wave function can be calculated 
efficiently as a determinant 0. As for the unpolarized 
system, we employ the fixed-node algorithm to avoid the 
fermion sign problem. This yields an upper bound to the 
energy for this system, the parameters in the BCS pair 
function <f)(rij) are (approximately) optimized to yield 
the best fixed-node energy. We include Jastrow correla- 
tions between anti-parallel and parallel spins in the trial 
wave function to minimize the statistical errors Q . These 
do not affect the energy, however, as they do not change 
the nodal surfaces where = 0. 

The quasiparticle spectrum at kpa = 00 is displayed 
in Figure □ The BCS prediction at k F a = 00 is shown 
for comparison. The QMC points are calculated by com- 
puting E k (N + 1) - [E (N) + E (N + 2)]/2 at constant 
density, where Ek(N + l) is the energy of the state of mo- 
mentum k with 1 unpaired and N paired particles. The 
Eq are the ground states of the N and N + 2 particle 
systems. 

Note that the minimum is at a momentum significantly 
less than the Fermi momentum. For larger coupling, the 
minimum in the dispersion will continue to trend towards 
lower momenta. This is apparent in the figure from the 
two sets of QMC results. From the QMC calculations we 
extract £ = 0.42(1) and a gap of 0.84(4)£ FG , or j3 = 1.2, 
somewhat smaller than earlier results. 

In order to understand if the BCS phase can support a 
finite polarization, it is important to also study the inter- 
action between the (polarized) quasi-particles. We have 
used QMC techniques to search for the ground state en- 
ergy as a function of polarization. The lowest variational 
energy in each case is found by filling the states at the 
minima of the quasiparticle spectra. 

The results are shown in Fig. |3| These results demon- 
strate that the quasi-particles at small polarization are 
nearly non-interacting. This would be expected at small 
polarizations since the pair size is expected to be of the 
order of the inter-particle separation. The solid points 
are the QMC calculations with finite polarization. The 
integers next to these solid points indicate the momen- 
tum shells n 2 , where k 2 = n 2 (2n/L) 2 = (n 2 x + n 2 + 
n 2 )(2ir/L) 2 , filled in the trial wave function. Calculations 
were performed for 66 particles at various polarizations. 

The open symbols represent the sum of single-particle 
energies obtained from the single-particle dispersion cal- 
culations of Figure|5] At small polarizations this is nearly 
degenerate with the full calculation. The solid line in 
the figure is the energy of a phase separated unpolarized 
paired phase and a fully polarized Fermi Gas. 

The calculations indicate that the mixed and homo- 
geneous phases are essentially degenerate at kpa = 00. 
Our calculations reproduce momentum distributions sim- 
ilar to those proposed for the gapless SC state, the states 
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Polarization 



that we have not analyzed other exotic possibilities such 
as the LOFF state which may be relevant. We are plan- 
ning on examining such possibilities in future calcula- 
tions, and also examining systems at large polarizations. 
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FIG. 3: Ground-state energy vs. polarization at fixed density 
and a = do. 



near the minimum of the dispersion are nearly fully occu- 
pied for one species, and nearly unoccupied for the other. 

For larger couplings (kpa > 0) at small polarizations, a 
transition to the homogeneous phase is expected to occur 
at some point from examining the behavior at very large 
couplings. Here the system can be thought of as a dilute 
mixed Fermi-Bose system where the bosons are tightly 
bound pairs of Fcrmions. The zero-temperature phase di- 
agram of such mixtures has been investigated by Viverit, 
Pethick, and Smith |23j . They find that at small fermion 
densities the system is homogeneous, while at larger den- 
sities it phase separates into either a pure fermion plus 
mixed or pure fermion plus pure boson phase. 

The results shown in Figure [3] are qualitatively simi- 
lar. Beyond a certain value of polarization, around 0.3 
in the figure, the homogeneous system is clearly higher 
in energy than a phase separated state consisting of an 
unpolarizcd BCS plus fully polarized fermions. In this 
regime it appears that the BCS component would con- 
tinue to support a smaller polarization. 

Conclusions: From our calculations it appears that at 
couplings near and beyond kpa = oo a homogeneous 
superfluid state with finite polarization and non-trivial 
gapless excitations may be accessible experimentally. We 
also find that at low polarization, the polarization is car- 
ried by quasi-particles that are nearly non-interacting 
and occupy momenta below kp. In cold atom experi- 
ments we expect the trap potential to favor the phase 
separated state, while finite-range effects, which are ex- 
pected to suppress the gap, and surface energy cost 
may help stabilize the homogeneous phase relative to the 
phase separated state. A naive interpretation of our re- 
sults would indicate that a gapless superfluid exists only 
at very strong coupling where A//i > 1 . However in more 
complex systems such as dense quark matter, charge neu- 
trality imposed by long-range forces plays an important 
role in this phase competition [L3J, |hj, |2J] . We also note 
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